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ABSTRACT 

Aims. To predict the cruise radiation environment related to future human missions to Mars, the correlation between solar modulation 
potential and the dose rate measured by the Radiation Assessment Detector (RAD) has been analyzed and empirical models have 
been employed to quantify this correlation. 

Methods. The instrument RAD, onboard Mars Science Laboratory’s (MSL) rover Curiosity, measures a broad spectrum of energetic 
particles along with the radiation dose rate during the 253-day cruise phase as well as on the surface of Mars. With these first ever 
measurements inside a spacecraft from Earth to Mars, RAD observed the impulsive enhancement of dose rate during solar particle 
events as well as a gradual evolution of the galactic cosmic ray (GCR) induced radiation dose rate due to the modulation of the primary 
GCR flux by the solar magnetic field, which correlates with long-term solar activities and heliospheric rotation. 

Results. We analyzed the dependence of the dose rate measured by RAD on solar modulation potentials and estimated the dose rate 
and dose equivalent under different solar modulation conditions. These estimations help us to have approximate predictions of the 
cruise radiation environment, such as the accumulated dose equivalent associated with future human missions to Mars. 

Conclusions. The predicted dose equivalent rate during solar maximum conditions could be as low as one-fourth of the current 
RAD cruise measurement. However, future measurements during solar maximum and minimum periods are essential to validate our 
estimations. 
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1. Introduction to MSL/RAD 


The Mars Science Laboratory (MSL) spacecraft, containing the 
Curiosity rover, was launched on November 26, 2011 and landed 
on Mars on August 6, 2012 after a 253-day cruise. T he Radia¬ 
tion Assessment Detector (RAD; iHassler et al.ll2012l) ) onboard 
Curiosity provided the first e ver radiation measu rements during 
its cruise from Earth t o Mars (Zeitlin et~^l2013h and now on the 
surface of the planet (IHassler et al. 2014t) bv measuring a broad 
spectrum of energetic particle radiation. The assessment of the 
radiation environment is fundamental for planning future human 
missions to Mars and evaluating the biologic effects of the radi¬ 
ation likely to be encountered on these kinds of missions. 


RAD has detected galactic cosmic rays (OCRs) and solar en¬ 
ergetic particles (SEPs), as well as secondary neutrons and other 
secondary particles created by the primary particles interact¬ 
ing with the spacecraft shielding. In a compact, low-mass, low- 
power instrument, RAD’s innovative design combines charged- 
and neutral-particle detection capabilities over a wide dynamic 
range. A schematic of the RAD sensor head (RSH) is shown 
in Eigure [H A deta iled overview of the detectors is given by 
IHassler et alJ (l2Q12l) . Eor charged particle detection, silicon de¬ 
tectors A and B are used in coincidence to define the acceptance 
angle of incoming particles. RAD measures mainly downward 
fluxes on the surface of Mars and it separates the charged par- 



Fig. 1. Schematic view (not to scale) of the RAD instrument, consist¬ 
ing of three silicon detectors (A, B, C, each having the thickness of 300 
pm), a caesium iodide scintillator (D with a height of 28 mm), and a 
plastic scintillator (E with a height of 18 mm). Both scintillators are sur¬ 
rounded by a plastic anticoincidence (El and E2). Eor detecting charged 
particles (blue). A, B, C, D, and E are used as a telescope. Neutral par¬ 
ticles (both neutrons and gammas, green) are detected in D and E using 
C, E l, and E2 as antico incidence. This figure is adapted from Eigure 
7 of iKbhler et al.l ( 120111) and originally from Eigure 3 of iPosner et all 
( 1200 ^ 
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tides in two categories: (a) particles stopping inside one of the 
detectors (B, C, D, or E), and (b) particles penetrating through 
the whole instrument and still creating energy depositions at the 
bottom detector of the sensor head. The total energy, charge, and 
mass of a stopping particle can be determined thus enabling its 
spectra to be reconstructed. The energy of a penetrating particle 
cannot, by definition, be easily recovered. However, the particle 
type and broad information about its energy can be retrieved us¬ 
ing the particle’s energy deposit pattern through the stack of de¬ 
tectors. The charged particle spectra m easured with RAD on th e 
surface of Gale crater are presented in lEhresmann et all (l2014l) . 
Eor neutral particle detection, both D and E detectors are, to dif¬ 
ferent degrees, sensitive to y-rays and neutrons. The separation 
and reconstruction of y-ray and neutron energy spectra on the 
Martian surface measu red by RAD have b een accomplished us¬ 
ing inversion methods (iKohler et al.ll2014li . 

2. Radiation environment and dose rate variations 
in transit to Mars 

In interplanetary space, OCRs contain 98% atomic nuclei and 
2% electrons, and the nuclei consist of about 87% p rotons, 12% 
helium, and only 1% heavier nuclei (ISimpsonlll983h . This com¬ 
position however varies slightly depending on the phase of the 
heliosphere modulation, which is related to both the 11-year so¬ 
lar activity cycle and the heliospheric rotatiorO]. Broadly speak¬ 
ing, the stronger the solar activity the higher would be the de¬ 
flection of interplanetary OCRs, resulting in a reduction of the 
GCR flux. This relationship can be often observed as an anti¬ 
correlation between neutron monitor counts (reflecting the GCR 
flux intensity at Earth) and the sunspot numbe r (a standard mea¬ 
surement of solar activity) dHeber et alj|2007h . This modulation 
is however energy-dependent: the flux of low-energy GCRs is 
more affected than that of the high-energy GCRs. 

On its journey to Mars, RAD started operating on Decem¬ 
ber 6 2011 and, with few interruptions, measured through July 
14 2012 when it was temporarily switched off in preparation for 
the landing. RAD collected more than seven months of first-time 
measurements of the radiation environment inside the shielding 
of the spacecraft from Earth to Mars, in principle, a similar en¬ 
vironment to future manned missions to the planet, thus pro- 
viding insight into the potential radiation hazard for astronauts 
(IZeitlin et al.ll2013]) . Because of the shielding of the spacecraft 
and internal structures, RAD measured a mix of primary and 
secondary particles. The latter are produced by primary particles 
via nuclear or electromagnetic interactions as they traverse the 
spacecraft. A simplified shielding model of the spacecraft devel¬ 
oped at JPL has been be used to calculate the shielding distri¬ 
bution as^een_by_RAJ^which is mounted to the top deck of the 
rover (IZeitlin et al.l2013h . Shielding around the RAD instrument 
during cruise was complex: most of the solid angle was lightly 
shielded with a column density smaller than 10 g/cm^, while the 
rest was broadly distributed over a range of depths up to about 
100 g/cm^. 

RAD measures dose in t wo detectors: the sil icon detector B 
and the plastic scintillator E dHassler et al.ll201^ . The latter has 
a composition similar to that of human tissue and is also more 
sensitive to neutrons than silicon detectors. The bottom panel in 
Eigure|2]shows the dose rate measurements taken by both detec¬ 
tors (black for E and gray for B) from 9 December 2011 to 14 
July 2012 during the cruise to Mars. Eor a given incident flux, 

’ The Sun has a latitude-dependent sidereal rotation period between 
25.4 days near the equator and 30-1- days at the poles. 


the dose rate in silicon is generally less than the dose rate in 
plastic because of the comparatively larger ionization potential 
of silicon. There were generally 44 dose rate measurements per 
day and we have calculated the average values of the dose rate 
in each day, shown in the figure as red and magenta dots for 
detector E and B, respectively. During quiet periods in solar ac¬ 
tivity, the average GCR dose rate measured by the silicon and 
plastic scintillators wa s about 332 + 23 p Gy/day and 461 +92 
pGy/day, respectively (IZeitlin et al.ll2013l) . In comparison with 
the measured dose rate in plastic, the dose r ate in silicon can b e 
converted using a conversion factor of 1.38 (IZeitlin et al.ll2013l) . 
which approximately relates energy loss per unit of path length 
(dE/djc) in silicon to linear energy transfer (LET) in water. Af¬ 
ter conversion, the dose rate in water, as measured in the silicon 
detector B, was 458 + 32 pGy/day, which is comparable to mea¬ 
surements in the plastic scintillator within uncertainties B 

The dose rate during the few SEP events (marked as the 
yellow-shaded areas) were as high as >10000 pGy/day. Eor a 
close visualization of the variations of the GCR-driven dose rate 
during the solar quiet periods, we omit the peak values of the 
SEPs in this figure. A zoomed-out figure containing the cruise 
dose rate meas urements, as well as the peak doses of the SEPs 
can be found in IZeitlin et al.l (l2013l) . 

2.1. Neutron monitor count rates and solar modulation 
strength 

Neutron monitor count rates probe the intensity of primary 
cosmic rays above certain thresholds in energy and magnetic 
rigidity, and therefore the measured flux depends on the lo¬ 
cation and sensitivity of the monitor. A commonly used pa- 
ra meter of heliospheric modulation is the modul ation potential 
<t> (iGleeson & Axfordlll968l: lUsoskin et al.ll2005l) . which corre¬ 
sponds to the mean energy loss of a cosmic ray particle inside 
the heliosphere and is often used to parametrize the modula¬ 
tion of the GCR spectrum. The parameter <1) is defined in a 
spherically symmetri c and steady-state case for the Earth’s orbit. 
lUsoskin et al.l (l2002l) used the following method to estimate <1) 
using the measured neutron monitor count rates: first, a stochas¬ 
tic simulation is used to calculate the galactic cosmic ray spectra 
at the Earth’s orbit for different yalues of <1); second, these spectra 
are conyoluted with the specific response function of the neutron 
monitor to generate the expected neutron monitor count rates 
giyen a certain yalue of <1); and Anally, the modulation strength 
O can be deriyed inyerting this forward-modeled relation from 
the actually recorded neutron monitor count rates. An estimated 
function of deriyation of <1) (in the unit of MV) from the mea¬ 
sured count rates (CR nm, in the unit of c ounts per minute) is 
giyen as the following (lUsoskin et al.ll20()^ : 


where A and B are modeled yalues that depend on the geograph¬ 
ical coordinates and hence the geomagnetic cutoff rigidity of the 
mo nitor stations . 

iPosner et al.l (l2013l) haye shown that during most of the 
cruise phase the MSL spacecraft was well connected with the 
Earth yia the Parker spiral magnetic held, a circumstance now 

^ Subsequent to IZeitlin et all ll2013li . we have re-evaluated the conver¬ 
sion factor used to obtain tissue dose from dose in silicon, and this yields 
the slightly lower dose rate (458 pGy/day) shown here compared to 
the earlier reported value (481 pGy/day). This is also the case for dose 
equivalent rate as shown later. 
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Day of Year 2012 


Fig. 2. Top: blue line shows the Oulu neutron monitor count rate 
recorded during the time of MSL’s cruise to Mars and the daily aver¬ 
aged values are shown as magenta dots. The hlack line shows the solar 
modulation potential O derived from the neutron monitor count rates 
and its daily average is shown as red dots. The RAD observed SEPs are 
highlighted in yellow shaded areas in which the data are not considered 
in the correlation study. Bottom: dose rate recorded hy RAD in silicon 
detector B (gray curve) and in plastic scintillator E (tissue-equivalent, 
hlack curve). Their daily averaged values are plotted as solid dots in 
magenta and red. 


referred to as the Hohmann-Parker effect. As a consequence, it 
was shown that the cosmic ray background at MSL and the Solar 
and Heliospheric Observatory (SOHO) near Earth were highly 
correlated. Therefore, for the evaluation of the solar modulation 
of GCR-induced MSL/RAD measurements, we can take approx¬ 
imately the modulation potential calibrated from measurements 
of the Oulu Neutron Monitor on Earth. The top panel of Eigure|2] 
shows both the count rates (blue, in the unit of counts per minute) 
recorded by the Oulu neutron monitor^ and the calibrated solar 
modulation potential O (black, in the unit of MV). We opted to 
use Oulu data with a resolution of three hours and for the current 
variation study, we calculated the daily average values of both 
neutron count rates and d) shown as solid dots in magenta and 
red, respectively. The daily variations of the measurements are 
folded into the sigma errors of the daily-average values. 

2.2. Correlation between solar modulation and cruise dose 
rate 

Impacts of solar events are often observed in the form of Eorbush 
decreases in the OCR flux due to the enhancement of t he inter¬ 
plane t ary magnetic field sweeping away more OCRs (iForbushI 
Il938t IZeitlin et al.ll2Q13l) . During SEPs and Eorbush decreases, 
when the GCR spectrum is distorted by the nonpotential en¬ 
hanced solar fields, the evaluations of <t) bear large uncertainties. 
Several Eorbush decreases and SEPs detected by the Oulu neu¬ 


^ The Oulu count rate data were obtained from 
http://cosmicrays.oulu.fi/ and the pressure effect was corrected. 


tron monitor and MSL/RAD are highlighted as yellow-shaded 
areas in Figure |2l They do not necessarily match precisely in 
time because of the spatial variability and time delay of solar 
events passing the two instruments. For the current correlation 
study of GCR-induced dose rate and solar modulations, we ex¬ 
clude the data collected during these five highlighted solar events 
since these measurements are due to both modulated GCRs and 
SEPs. During the quiet period, we calculated the daily averaged 
neutron count rate CRnm, modulation potential <1), as well as the 
daily averaged RAD dose rate shown in the bottom panel of Fig- 
urel^as red and magenta dots for E and B detector, respectively. 
The sigma error of each daily average value is taken to be the 
standard deviation of all the measurements within that day. The 
daily average data clearly indicate an anticorrelation between the 
variation of the RAD dose rate and the solar modulation poten¬ 
tial d). The correlation coefficient between two variables can be 
obtained as the covariance of the two variables divided by the 
product of their standard deviations, giving a value between -i-1 
and -1 inclusive, where 1 is total positive correlati on, 0 is no cor - 
relation, and -1 is total negative correlation (e.g.. lStigleijll989l) . 
We obtained the correlation coefficient between the daily aver¬ 
aged solar modulation O and the RAD dose rate (during periods 
with no SEPs) to be -0.80 for detector B and -0.77 for detector 
E, indicating a very good negative correlation between the solar 
modulation and the measured GCR induced dose rate. 

The radial gradient of galac tic cosmic ray s in the inner he¬ 
liosphere has been studied by iGieseler et al.l (l2008h and they 
estimated GCR flux increases at a radial gradient rate of 4.5% 
per AU in the inner heliosphere. Considering that the spacecraft 
moved progressively outward in the solar system from 1 AU to 
about 1.5 AU, we corrected this gradient increase on the mea¬ 
sured dose rate. This effect is however very small (within 2% of 
correction) and has only affected slightly the analysis that fol¬ 
lows. 

A quantitative study of the anticorrelation is challenging be¬ 
cause of the lack of analytic models. We employed two separate 
empirical models to describe this anticorrelation. 

1. A simple linear correlation described by a first order polyno¬ 
mial function is used to fit the daily average data, i.e.. 


D =/?o 


( 2 ) 


where D is the RAD measured dose rate (either by de¬ 
tector B or E) in the unit of [xGy/day and Pi in the unit 
of pGy/day/MV represents the linear relation between the 
changing of the solar modulation and the variance of ob¬ 
served dose rate. This linear correlation can be easily fit¬ 
ted by a simple regression method, which fits a straight line 
through the set of measurements minimizing the sum of 
squared residuals (i.e., the chi square function) of the model. 
Sigma errors of the data can also be included in the fitting as 
a weighting factor. The fitted parameters and their standard 
deviations are j3o - 626+48 [.iGy/day and /3i - -0.39 ± 0.07 
pGy/day/MV for dose rate in B and /So = 758 +27 qGy/day 
and jSi - -0.44 + 0.04 [xGy/day/MV for dose rate in E. Given 
the small error bars of the fitted parameters (< 15%), a linear 
function seems to be a good approximation for the range of 
the measured modulation, i.e., between 550 MV and 800 
MV. 

2. We employ an empirical model, whic h is often used in th e 
analysis of neutron count rates (e.g.. lUsoskin et ani201 Ih . 
i.e.. 


D 


ao + 


O'! 

<h + a2 ’ 


( 3 ) 
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Fig. 3. Data and fittings processed through bootstrap Monte Carlo sim¬ 
ulations of the variations of the RAD dose rate and the solar modulation 
O during the cruise to Mars.The red/magenta lines represent the linear 
fits (Eq. |2j of RAD measured dose rate in detector E/B versus solar 
modulation potential. The green/cyan lines show the results of the non¬ 
linear fits (Eq.[3)- 


with the requirement of a 2 > 0. Otherwise, (t>+a 2 - 0 would 
lead to nonphysical solution of the problem. Also by mini¬ 
mizing the chi square function, we fitted the above function 
and obtained the parameters as ao - 94.4 + 2.6 qGy/day, 
ai = (1.8 + 0.02 )x 105MV ■ (xGy/day and a 2 = 0.0002 + 1.0 
MV for dose rate in B and ao = 180.6 + 0.7 [xGy/day, 
ai = (1.9 + 0.005) X lO^MV ■ (xGy/day and a 2 = 0.0008 + 1.4 
MV for dose rate in E. The error bars for both ao and ai are 
below 3% while the relative error for a 2 is much bigger be¬ 
cause of the small absolute values of a 2 , which indicates that 
the parameter a 2 could be probably ignored in the model fit¬ 
ted by the current measurement. 

2.3. Estimates of the cruise radiation environment under 
different soiar moduiation conditions 

It is difficult to derive analytically whether the above models 
would serve as good approximations when solar modulation <D> 
is outside the range between 550 and 800 MV of the studied 
case. In general, solar modulation affects particles with different 
energies in different degrees, i.e., <1) has a str onger/weaker effec t 
on GCR particles with lower/higher energies dHeber et al.l2007h . 
For any given species, a lower-energy ion contributes more to 
dose than does a higher-energy ion; we therefore expect that so¬ 
lar modulation has a more pronounced effect on dose than it 
does on integral fluxes. Furthermore, the RAD dose rate dur¬ 
ing cruise was measured under the spacecraft shielding, which is 
more effective against lower energy particles that are more influ¬ 
enced by the variation of fl). A more detailed quantitative analy¬ 
sis could be possibly car ried out in the future u sing a GEANT4 
Monte Carlo simulation dAgostinelli et al.ir2003h considering the 
energy-dependent solar modulations and the energy-dependent 
detector responses as well as the complex shielding model of the 
spacecraft. Nevertheless, based on our empirical regression fit¬ 
tings of the above two models, we could extrapolate our results 
to a wider range and make rough estimations of the dose rate 
values under different solar modulation conditions. 



Fig. 4. Annual average of reconstructed modulation potential O since 
1937 (gray line and dots with units on the left axis) and estimated annual 
dose equivalent rate derived from different models (with units on the 
right axis). Cyan/Green line shows the dose equivalent values derived 
from the nonlinear model (Eq.[^ based on detector B/E measurements. 
Magenta/Red line represents the dose equivalent values derived from 
the linear model (Eq.|2} based on detector B/E measurements. Black, 
yellow, and blue dots represent the data within the range of modula¬ 
tion fl) between 550 and 810 MV which the RAD cruise measurements 
covered. 

A simple error propagation method based on the fitted pa¬ 
rameters and their errors is however not reliable when extrapo¬ 
lating the function away from the measured data. To derive the 
robust estimates of the uncertainties due to both measurements 
and extrap olations, we h ave carried out a bootstrap Monte Carlo 
simulation (lEfronll98lh . Simulated data sets are generated using 
the uncertainty range of the measured data, which are then used 
for fitting the model (according to either Eq.|2]or Eq.O and ex¬ 
trapolating to a wider range of <D> to obtain an accurate estimate 
of the uncertainties in the extrapolated data. 

The regression fittings of the above two models and the ex¬ 
trapolations of the fittings (over a range of <6 from 250 to 1200 
MV), using bootstrap Monte Carlo simulations, are shown in 
Figure |3] As expected, the uncertainty of the extrapolated dose 
rate increases when the extrapolation is further away from the 
actual measurements. The differences of dose rate between two 
models are bigger for extrapolated points, especially for small 
solar modulation, i.e., solar minimum case. For instance, at 
O = 25QMV , the discrepancy between the two models is as large 
as 300 pGy/day for the plastic detector E. Nonetheless, prefer¬ 
ence of a better model is difficult to reach because of the limited 
statistics and range of the current data set. Future measurements 
over a wider range of solar modulations are necessary prior to 
human flights to Mars to help improve and determine the mod¬ 
els. 

2.4. Predictions of the dose equivaient rate 

The Oulu neutron monitor (lUsoskin et al.ll20flli published the 
annual values of reconstructed modulation potential <1) since 
1937 and, accordingly, we are able to calculate the extrapo¬ 
lated dose rate values. We also correlated the monthly <I) values 
and monthly sunspot number provided by the Royal Observa¬ 
tory of Belgium, Brussels (http;//sidc.oma.be/silso/dataflles) us- 
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ing the same method as described in the last section. Based on 
the sunspot predictions in the next five years (2015 to 2019) pre¬ 
pared by the U.S. Dept, of Commerce, NOAA, Space Weather 
Prediction Center (SWPC), we estimate the corresponding mod¬ 
ulation potential as shown as the last five points in Figure |4] 

Dose equivalent, in Sieverts (Sv), defined by the Interna¬ 
tional Commission on Radiological Protection, is taken to be 
proportional to the risk of lifetime ca ncer induction via popu¬ 
lations studies (e.g. JZeitlin et al.ll2013h . It can be approximately 
estimated to be <Q> x D, where <Q> is the average quality 
factor, which is a conventional parameter for radiation risk es¬ 
timation. A value of 3.82 + 0.2 5 was measured during MSL’s 
cruise phase dZeitlin et alJl2013]l . Multiplying the measured tis¬ 
sue dose rate with the average quality factor yields the GCR 
dose equivalent rate. For the dose rate measured by detector B, a 
silicon-water conversion factor of 1.38 has to be applied first as 
described in Section|3 

The extrapolated dose rate estimations at different modula¬ 
tion potential values can therefore be used to evaluate the dose 
equivalent rate as shown in Figure |4] The annual average of re¬ 
constructed modulation potential (D shows a clear 11-year cycle 
following the solar cycle and it varies over a wide range from 
less than 250 MV to more than 1200 MV; the estimated annual 
dose equivalent rate ranges from about 0.5 mSv/day to about 5 
mSv/day and has a clear anticorrelation with O, as expected from 
both of our models (Eqs.l^andO; stronger solar modulations 
lead to a decrease of dose equivalent rate and, at O > 1000 MV, 
the dose equivalent rate can be as low as < 1 mSv /day within 
the uncertainti es, approximately half of the average RAD cruise 
measurement dZeitlin et al.ll2013h : 1.75 + 0.30 mSv/day; at solar 
maximum and minimum when the modulation potential is fur¬ 
ther away from the measured range, the discrepancy between the 
two models increases; the linear model often predicts a smaller 
dose rate than the nonlinear model as also shown in Figure [3j 
large uncertainties emerge at low solar modulation potentials es¬ 
pecially for the nonlinear model and future measurements over 
solar minimum periods prior to human flights to Mars are essen¬ 
tial for improving the predictions at low modulation potentials. 


3. Discussions and conciusions 

We presented the dose rate data collected by MSL/RAD during 
its cruise phase from Earth to Mars and analyzed its variations 
during solar quiet periods. The variation of the GCR-induced 
dose rate at MSL occurs nearly simultaneously with the varia¬ 
tion in neutron monitor count rate at Earth and was mainly driven 
by changes in heliospheric conditions. The simultaneity of cos¬ 
mic ray count rate changes at MSL and Earth is attributed to 
Hohm ann-Parker effect conditions as described in iPosner et al.l 
(l2013h . We show here that dose rate at MSL during the Mars 
transfer can be described by a solar modulation factor (<t)) calcu¬ 
lated from the neutron monitor count rate on Earth. The results 
show a clear anticorrelation trend, with a correlation coefficient 
of » -0.80 between <I) and the measured dose rate, indicating 
that when the solar activity is stronger GCR particles are attenu¬ 
ated more resulting in a decrease in the dose rate measurements. 
The measurement from DOSIS onboard ISS has also shown 
an anticorrelation between solar modulation and the dose rate 
in the silicon detector (through personal communications with 
Sonke Burmeister). However, the results are not directly compa¬ 
rable with RAD’s measurement because ISS is in low Earth orbit 
(LEO) where Earth itself and its magnetic field provide shielding 
against GCRs. 


A quantitative study of this anticorrelation has been carried 
out employing two empirical models, which fit the data equally 
well as shown by the small standard deviation of the fitted line. 
We further extrapolated the regression fittings to a wider range 
of solar modulations to make approximate empirical predictions 
of dose rate and dose equivalent rate under different modulation 
conditions. We employed a bootstrap Monte Carlo simulation 
to include possible propagated errors due to the uncertainties of 
both measurements and extrapolations. The final derived dose 
equivalent rate at various modulation potential values are shown 
with a clear anticorrelation. The predicted dose equivalent rate 
during solar maximum conditions, where <I) ~ 1200 MV was 
found to be as low as » 0.5-1.5 mSv/day (as within the uncer¬ 
tainties of both models), which is considerably lower than the 
RAD cruise measurement 1.75+0.30 mSv/day. The discrepancy 
of the two fitted models, as well as the extrapolation uncertain¬ 
ties, are significant at low solar modulation potentials and future 
measurements during solar minimum periods are necessary for 
improving the predictions at this range. 

We also estimated the solar modulation potential in the next 
five years, using its correlation with sunspot numbers and em¬ 
ploying the same method of correlating and extrapolating dose 
rate and O. As the predicted solar modulation will decline in the 
future years, we predict a trend of increasing dose equivalent rate 
(between 1.7 and 3 mSv/day) from year 2015 until 2020. 

Total mission GCR dose equivalent can be estimated based 
on our measurements and predictions. Considering a similar 
shielding condition, assuming a 180-day, one-way duration as a 
typica l NASA’s “Design Reference” Mars mission (iDrake et al.l 
1201 Oh . we could estimate a crew taking a 360-day round trip to 
receive about 360 + 180 mSv (the error bar combines the sigma 
from both models) from GCRs under a high solar modulation 
condition (O s; 1200 MV). The fastest round trip with on-orbit 
staging and existing propulsion technologies has been estimated 
to be a 195-day trip (120 days out, 75 days back with an extra , 
e.g., 14 days on the surface), as described bv iFolta et al.l (l2012h . 
This would result in an even smaller GCR-induced cruise dose 
equivalent during solar maximum; 195 + 98 mSv. 

However, this fast round-trip option requires extra propul¬ 
sion, which may result in a reduction of payload mass and a 
change of shielding conditions. It is difficult to estimate how 
reduced shielding would affect the dose rate, but it is probable 
that particles would traverse through more easily, heavy ions 
less likely to break up into lighter ions, and the production of 
secondary particles, such as neutrons, reduced. The dose rate 
and the quality factor are expected to be modified, both fold¬ 
ing into the changes of dose equivalent. Monte Carlo GEANT4 
simulations could be carried out for investigating this compli¬ 
cated process where realistic shielding geometry and materials 
can be included. In the case of the MSL spacecraf t, the shield¬ 
ing co ndition around RAD was highly nonuniform dZeitlin et al.l 
l2013h ; nearly 30% of incoming particle trajectories within the 
field of view of RAD from above (coincidence of detector A 
and B in Eigure[TJ were only lightly shielded (< 5 g cm“^), the 
50% fraction from below were shielded by 8-10 g cm“^, and 
the remaining 20% were heavily shielded (> 20 g cm“^), giv¬ 
ing an average depth of roughly 16 g cm“^. It is highly proba¬ 
ble that the crewed mission would not benefit from a perfectly 
uniform shielding. Also, the shielding materials for MSL have 
been rather complex; the spacecraft has a mix of materials that 
included hydrazine, mylar (parachute material), as well as alu¬ 
minum and carbon composite (spacecraft shell). 

However, we obtained the correlation between solar modu¬ 
lation potential <t) and dose rate based on a very good magnetic 
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connection between Earth and the spacecraft, allowing the di¬ 
rect application of the Neutron Monitor data collected at Earth. 
In the case of a weak magnetic connection between Earth and 
the spacecraft such as in a rapid Mars transit scenario, the local 
solar modulation potential <t> at the spacecraft location should 
be considered. Our current analysis simply assumes that the so¬ 
lar modulation measured at Earth is representative for that at the 
Earth orbit and we include the radial gradient of this modulation. 
Euture work may also consider the longitudinal difference of the 
modulation across the Parker spirals. 

It should also to be stressed that additional contributions 
of dose rate and dose equivalent rate by SEPs can differ sig¬ 
nificantly from our current measurements since the frequency 
and intensity of SEP events are highly variable. The MSL/RAD 
data were obtained under the rising and maximum phase of 
an unusually weak solar cycle. Only five SEP events were ob- 
ser ved during the 253 -day cruise to Mars as already reported 
by IZeitlin et alj (l2013h : two from 23 to 29 January with total 
dose equivalent being 4.0 mSv, two from 7 to 15 March with 
total dose equivalent being 19.5 mSv, and one on 17 May of 
2012 with total dose equivalent being 1.2 mSv. The summed 
dose equivalent of the five observed SEPs is 24.7 mSv, roughly 
5% of the contribution to the total cruise dose equivalent (» 466 
mSv). During the peak of a solar maximum, where the solar ac¬ 
tivity is much more intensive, one would generally expect a big¬ 
ger contribution of SEP-related dose equivalent. Nonetheless, a 
shorter trip (360-day or 195-day return) would also reduce the 
risk of encountering many SEP events. If we assume that on a 
180-day trip during solar maximum peak a crew would receive 
four times the total SEP-contributed dose equivalent, i.e., 100 

mSv, a round trip to Mars would result in a total dose equiva¬ 
lent of 560 + 180 mSv, still smaller than the estimation based 
on the current solar modulation conditio ns, i.e., 662 ± 108 mSv 
dZeitlin et al.ll2013t iHassler et alJl2014h . A 195-day round trip 
would of course further reduce our dose level estimations. These 
estimations are slightly less than the safe upper limit for 30- to 
60-years old, nonsmoking females (600-1000 mSv) and males 
(800-1200 mSv) given by th e Central estimates of dose limits 
dCucinotta & Chapi^lll201 lb . However, given the uncertainties 
in the models, extra caution would be necessary in planning a 
human mission during solar maximum. Thus, future measure¬ 
ments over maximum solar modulation potentials are essential 
to refine our predictions of the reduced dose rate at this range. 
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